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Recent polarized neutron scattering experiments on Y Ba^CuzOi have revealed a sharp spectral 
peak at the (tt, vr) in reciprocal lattice centered around the energy transfer of 41 meV . We offer 
a theoretical explanation of this remarkable experiment in terms of a new collective mode in the 
particle particle channel of the Hubbard model. This collective mode yields valuable information 
about the symmetry of the superconducting gap. 
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in 



> 
o 
vo 
o 

o 
in 

o> 

S 

i 

C 

O 

o 



Recently, both unpolarizcd and polarized neutron 
scattering experiments have been performed on the 
Y Ba-iCu^O'j high T c superconductors |^-|^]. In particu- 
lar, polarized neutron experiment shows an extremely 
sharp spectral feature in the spin flip channel. This fea- 
ture is centered around (jr, 7r) in reciprocal space, and 
peaked at 41 meV with a width narrower than the in- 
strumental resolution. This feature also has an inter- 
esting temperature dependence. In the experiment by 
Mook et al, while it exists above the superconducting 
transition temperature of T c — 92AK , its intensity scales 
like the superfluid density below the transition. More re- 
cently, Fong et al performed detailed spin unpolarizcd 
neutron experiments with a careful subtraction of the 
phonon background. They found that the 41 meV mode 
disappears above the superconducting transition temper- 
ature. 

In this letter, we offer a theoretical explanation of this 
remarkable experiment. We first show that for a general 
class of tight binding Hamiltonian including the Hubbard 
and the t — J model, there exist well-defined collective 
modes in the particle particle channel centered around 
momentum (w, tt). The spin quantum number of this ex- 
citation can either be a singlet or a triplet. The singlet 
excitation has been discussed by Yang Q and one of us 
P|, and is in fact an exact eigenstate of the Hubbard 
model. Normally, collective excitations in the particle 
particle channel are inaccessible experimentally. How- 
ever, one of us || argued that if the ground state of the 
model in consideration is superconducting, one can cou- 
ple to it through a particle hole excitation, because the 
BCS condensate is a coherent mixture of particles and 
holes. Based on this consideration, one of us |fj predicted 
a possible new collective mode of the high T c supercon- 
ductors. It is a spin singlet excitation peaked at (7r,7r), 
has a well defined energy of [/ — 2/i and its intensity scales 
like the superfluid density. Possibly because it is hard to 
distinguish it from other excitations in the system, this 
mode has not yet been detected experimentally. 

However, the basic arguments can be easily generalized 
from the singlet to the triplet case. Besides the above 
mentioned collective mode in the singlet particle particle 



channel, there also exists a well defined collective mode 
in the triplet channel near total momentum (7r,7r). This 
is true for a large class of tight binding models, like the 
Hubbard or the t — J model. The energy spectrum of 
a non-interacting pair of particles or holes generally con- 
sists of a continuum labeled by their relative momentum. 
However, for tight binding models, this continuum col- 
lapses to a point when the total momentum of the pair 
is (it, it). (See Figure 1). 
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FIG. 1. The energy spectrum along the (ir,n) direction. 
The dots correspond to the antibonding state and the solid 
line is the edge of the continuous spectrum. Here the numer- 
ical calculations were done for J = t and n = 0.85. As 8 we 
denoted ^j. 

A triplet pair generally has a repulsive interaction 
when placed on the neighboring site. Because of the 
collapse of the particle particle continuum, this repul- 
sive interaction leads to an anti-bound state near total 
momentum (7r, 7r, ...) in any space dimensions. (See Fig- 
ure 1). This anti-bound triplet state manifests itself as 
a collective excitation of the many-body system, or as a 
pole in the Greens function of the particle particle chan- 
nel. Most physical probes do not couple to this channel. 
However, if the model in consideration is superconduct- 
ing (which is an assumption of our theory), then a spin 
flip scattering of the neutron can couple directly to the 
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triplet particle particle excitation. Physically this pro- 
cess is nothing but the spin flip scattering of a Cooper 
pair in the RCS condensate. There is an overwhelming 
evidence that the Cooper pair in the high T c materials is 
a spin singlet with total momentum zero Therefore 
a spin flip scattering of the Cooper pair creates a triplet 
state. The orbital angular momentum or the parity se- 
lection rule forbids such a transition if the momentum 
loss of the neutron is zero, on the other hand, it can be 
shown simply that the matrix element of this coupling is 
maximal when the momentum loss of the neutron is at 

(7T,7r). 

Based on the above reasoning, we shall interpret the 
sharp spectral feature observed in the polarized neutron 
experiment in terms of the triplet collective mode in the 
particle particle channel. We begin with some analytic 
calculations of the properties of this collective mode and 
its contribution to the dynamical spin spin correlation 
function. Subsequently, we compare our predictions with 
a number of characteristic features of the resonance ob- 
served in experiments. Finally we point out some difficul- 
ties with our interpretation and their possible resolution. 

Our anti-bound state is different from the excitons of a 
superconductor considered by Bardasis and Schrieffer . 
These excitons form because of an attractive potential in 
a given angular momentum channel, and they exists in- 
side the superconducting gap and near total momentum 
q = 0. Our model is also different from the magnetic sus- 
ceptibility in the superconducting state computed using 
RPA H , since it involves multiple scattering in the par- 
ticle hole channel. More recently, Bulut and Scalapino 
considered a model of the bi-layer superconductors 
and argued that a dynamic nesting from the bonding to 
anti-bonding Fermi surface could give rise to a collective 
resonance at (it, tt). Our main difference lies in the fact 
that their resonance is in the particle hole channel and 
as such it can exist well above T c . The bi-layer band 
structure plays a crucial role in their model while it is 
irrelevant in our case. 

We consider the following model defined on a two di- 
mensional square lattice: 

H = - t c\a c 3<y + JJ2^- Sj 

<ij> <ij> 

+ U ^2 n^nn - At J] c ^ ClCT W 

i i 

In the limit U — > oo, we recover the t — J model. On the 
other hand, if we keep U small, we can regard the above 
model as an effective Hamiltonian of the weak coupling 
Hubbard model, where the J term arises from the param- 
agnon interaction on the nearest neighbor sites. To study 
the spectrum in the triplet particle particle channel, we 
consider the operator 

Oj = £/<r(P)4+ 8 T c -PT ( 2 ) 
p 

which describes a pair of particles with center of mass 
momentum q and relative momentum p. The equation of 



motion for this operator is given by its commutator with 
the Hamiltonian, [H, Oj] . If we factorize the resulting 
commutator in terms of Oj and the expectation values 
of the density riko- , we obtain: 

[H,0\] = 

J2(t P +q + e -p)/ g b)Cp+?T c -pT + 2 ( Un l ~ P)°l 
p 

p p' 

x iv(p-p') - v(p + p' + <?)) (! ~ ?V-«T _ n p't) 

-|rE/^) c P+ 9 T ct -Pt 
v 

x F £y u p' fa(P + P') + Vip' -p-d)) (3) 

where n a = -k n& CT is the average density of electrons 
with spin er, N is the total number of lattice sites and and 
7](p) = X)s ex P[*P • o] is a geometrical factor coming from 
the summation over the nearest neighbors. In this equa- 
tion, the first term describes the kinetic energy of the pair 
of particles in consideration. For tight binding models 
with nearest neighbor hopping, efc = — 2t(cos k x + cos k y ). 
In this case, the kinetic energy of the pair vanishes when 
the total momentum q — (tt,tt). The second term de- 
scribes the Hartree interaction of the spin up pair with 
the average density of the down spins in the background 
and the chemical potential energy of adding a pair of 
particles. In the large U limit, the leading contribution 
to the chemical potential is given by Un/2. In fact, for 
Hubbard model at half filling, i.e. at n = 1, /i = U/2 
is an exact relation. Therefore, the second terms can- 
cels in the leading order in U, and reaches a finite limit 
as U — > oo. The third term gives the multiple scatter- 
ing of the particles with each other on a restricted phase 
space due to the filled Fermi sea. The last term describes 
the Fock self-energy of the quasi-particles. The method 
of factorizing the operator equation of motion is equiva- 
lent to the T matrix approximation in the diagrammatic 
calculations, it is exact in the low density limit. The col- 
lective mode for the triplet particle particle excitation is 
obtained by equating the right hand side of (3) to E q O^, 
the resulting eigenvalue equation is given by 

J sin 2 (p Q )(l-ng pT -ng +pT ) 

1 = > — (4) 

2JV^E«+4i~[cos(p x )cos(f ) +C os(Ps,)cos(f )] 

with a being x or y, i = t + ^2 p n P i cos(p x ) and 

E q = E q + 2(Uni — fi). It is straightforward to see that 
a discrete eigenvalue exists at q = (tt, 7r), with energy of 

E=^(l-j^J2n p cos 2 Px ) (5) 
p 

Discrete eigenvalues also exist for a finite range in q space. 
There are in general two eigenvalues for each q, one cor- 
responds to a p wave pair oriented in the x direction, 
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while the other is oriented in the y direction. The two 
eigenvalues are degenerate only along the diagonal when 
Qx = Qy The above eigenvalue equation can easily be 
solved numerically and the dispersion of the collective 
mode is shown in figure 1 . £ is basically the renormalized 
hopping matrix element of the quasi-particle. The range 
in momentum space over which the collective mode exists 
depends on the ratio of J ft. Numerous numerical calcula- 
tions indicate that the quasi-particle bandwidth is of the 
order of J JlOj ] . Here we take a semi-phenomenological 
approach and choose t rather than t as a free parameter. 
For a ratio of J/t = 1, and n — 0.85 we see that the 
collective mode exists over of the momentum space. 

In the experiments with the polarized neutrons scat- 
tering one measures the dynamic spin-spin correlation 
function 



where A q is the normalization factor of the wave function 



S{q,w) = Y,\^\S\\0) \ 2 S(co-u; n0 ) 



(0) 



where | > and | n > are the ground and excited states 
of the system and 5| = J2 P c P + q -\ c pl- Using the operator 
equation [H, O qQ \ = E q Oj a we can construct a class of 
approximate excited states of the Hubbard Hamiltonian 
as | n) = O qa | 0). The same operator equation for O g a 
shows that O q _ a | 0) = 0. With these relations, we can 
calculate the contribution of this approximate eigenstate 
to S(q, to) at zero temperature: 



= £<0|O 9 , Q S, 



0) | 2 5(u-EZ) + Y, 

a. n' 

£<0| [O q ^S\] |0) \ 2 5{u-E« 



E 



(7) 



where , denotes the contribution from states other 
than O q a | 0), and we used the fact that O q _ a 0) = 
to replace the product of two operators by their commu- 
tator. Evaluating the commutator we obtain 

So(q,u) = 

I (0 I 2^/^(p- hc_ p+it c p _ u | 0) | 2 5{u>-E«) 



4 E I E ti a) (p I 2 *(«> - K) 



(8) 



Summation over a corresponds to the two kinds of an- 
tibonding solutions. The overlap matrix element is fi- 
nite for both extended s wave and d wave pairing in the 
ground state. If we take the d wave symmetry of the order 
parameter A^. = Ao(cos(fc 2 ,) — cos(k y )) and, for example, 

(x) 

f q - a p wave oriented in the x direction. Then, 
<3\ 



p 
1 

X 



9 V % + 4 *~( C0S (P*) cos(f ) + cos( Py ) cos(f-)) 
2A . q x 



A 2 - 



E 



sm 2 (p x )(l - nq p - nq +p f 



7 \E q + 4t(cos{p x ) cos(^) + cos(p y ) cos(f-)) 



Finally we have 



S(q,u,T) = l6 



^[sin 2 (|)^-^) 



(9) 



+ sin 2 (|)^-^)](l + n B (^)) (10) 

where the last factor takes care of the finite tempera- 
ture in spin-spin correlation functions. One can easily 
see though, that the main temperature dependence of 
S(q,u>,T) comes from A (T). We also see that the in- 
tensity vanishes at zero external momentum and is max- 
imal at (tt, it). This fact can simply be attributed to 
the conservation of orbital angular momentum. Since 
the Cooper pair in the ground state has angular momen- 
tum either I = or I = 2, a finite momentum transfer 
is required for it to be excited to a state with angular 
momentum 1 = 1. The dispersion of two states is not 
very big giving rise to only a small splitting of two levels 
which never exceeds 15%. 

We also see that the spectral intensity is simply pro- 
portional to the BCS order parameter or the superfluid 
density. This is so because the BCS order parameter pro- 
vides the coupling from the particle hole channel to the 
particle particle channel This extra spectral weight 
at energy E q ^ a is transfered from the low energy sector, 
since the singlet BCS pairing removes the low energy spin 
fluctuation. Our theory is consistent with the fact that 
the 4lmeV peak intensity seems to scale with the su- 
perfluid density below T c . This conclusion agrees with 
the experimental results obtained by Fong et al but 
it is inconsistent with the fact that this mode does not 
vanish above T c in the experiment by Mook et al . At 
this moment, it is an experimental question that has to 
be resolved. We also note that our calculation is not 
fully self-consistent, and assumes that the mode energy 
does not change significantly when the system becomes 
superconducting. This is true for a variety of collective 
modes like the plasmon and the singlet particle particle 
mode j|,H (the so called 77 pair). But we do not have a 
general proof in the present case. A fully self-consistent 
calculation will be carried out in the future [ll| ] , but we 
anticipate that the basic physics remain unchanged. 

All the above discussions were restricted to the two 
dimensional CuO plane, which we model by the two di- 
mensional Hubbard model. However, it can be simply 
generalized to three dimensions. If one takes a three di- 
mensional Hubbard model, one finds that collective mode 
exists near (7r,7r, it) rather than (77, tt, 0). This is consis- 
tent with the experiment where the third component of 
the momentum transfer is also tt. 

The fact that the particle particle collective mode al- 
ways exists at (tt,tt, ..) is a special property of the tight 
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binding model on a bi-partite lattice. How would a next- 
nearest neighbor hopping term change the results? In 
this case, the dispersion relation is given by 



-2 1 (cos k x + cos k y ) — 2t ' (cos ( k x + k y ) + cos ( k x — k y ) ) 



where t' denotes the amplitude of the next-nearest neigh- 
bor hopping. The x and the y collective modes are mixed 
into the bonding sin(p K ) + sin(jp y ) and the antibonding 
sin(p x ) — shx(py) combination, and we obtain a matrix 
equation for the eigenvalues of these collective modes: 



det 



where z r . 



J 1 y-v sm(p a )sm(pp)z p 



2 N 



<WII=o (11) 



1 



nz +p and Sl pq 



<-pq — * '"j-p 
It is easy to see that the antibonding state E~ v always 
exists, while there is a critical coupling 



t' = 1.2 x 10" 4 J. 



(12) 



for the bonding state, so that for t' > t' cr , it ceases to ex- 
ist at (tt, tt). In real experiments t' and J are of the same 
order of magnitude, therefore, one can safely conclude 
that the anti-bonding state disappears into the contin- 
uum. This conclusion leads directly to a measurement 
of the symmetry of the superconducting gap. Below the 
superconducting transition temperature, the intensity of 
the anti-bonding state as measured in the neutron scat- 
tering experiment is proportional to 



E4 a) (p-5v 



?/2 



(13) 



where fq (p- 



oc sin(p x ) — sm(p y ). From this equation 



we sec immediately that the intensity is non vanishing if 
and only if the gap symmetry is of the d wave type. We 
therefore argue that the existence of the neutron reso- 
nance determines the symmetry of the pairing gap of the 
high T c superconductors to be of the d wave type, con- 
sistent with the theories wherepairing interaction arises 
from the spin fluctuations |12] -Q . It could also be consis- 
tent with more exotic possibilities of d x i_ y i +id xy pairing 
symmetry |Ie}| . 

We conclude that the basic features of the observed po- 
larized neutron scattering experiment can be explained 
in terms of a new particle particle collective mode in the 
Hubbard model. The energy of the mode at (tt, tt) is 
given by (|J) in the case of t' = and implicitly by ( |i"T|) for 
finite t' . Unfortunately, the energy of the mode depends 
on all three parameters /i, J and tf . While these param- 
eters can be fitted to the experimental value, the fitting 
itself is not a confirmation of our theory. The mode is 
centered around (tt, tt, tt) for two different reasons, both 
because the particle particle continuum at this momen- 
tum is minimal so that the anti-bound state could exist, 
and because of the conservation of angular momentum 
for exciting a singlet Cooper pair to a triplet state. This 
is exactly the momentum transfer of the excitation ob- 
served in experiment. The intensity of the mode scales 



with the superfluid density because the BCS pairing am- 
plitude is involved in converting a particle hole pair into 
a particle particle pair. Above T c , our simple calculation 
shows that the mode disappears, a much more sophisti- 
cated calculation is required to understand if the fluctu- 
ation pairing effects could give rise to a resonance even 
above T c . When the effects of next- nearest- neighbor hop- 
ping is included, only the anti-bonding sin(p x ) — sm(p y ) 
mode can exist and the region in momentum space for the 
existence of this collective mode becomes much narrower. 
However, the sharp collective mode might still exist as a 
broad resonance outside of this region. The antibond- 
ing collective mode only has an overlap with the d wave 
order parameter, and we conclude that the experimen- 
tal observation of the collective resonance in the neutron 
scattering experiment can only be consistent with the d 
wave symmetry of the pairing gap. 
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